Evidence of spontaneous spin polarization in the two-dimensional electron gas 
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Density-functional calculations using an exact exchange potential for a two-dimensional electron 
gas (2DEG) formed in a GaAs single quantum well predict the existence of a spin-polarized phase, 
when an excited subband becomes slightly populated. Direct experimental evidence is obtained from 
low temperature and low excitation-power photoluminescence (PL) spectra which display the sudden 
appearance of a sharp emission peak below the energy of the optical transition from the first excited 
electron subband upon its occupation. The behavior of this PL feature in magnetic fields applied 
in-plane as well as perpendicular to the 2DEG indicate the formation of spin-polarized domains 
in the excited subband with in-plane magnetization. For it speaks also the strong enhancement of 
exchange-vertex corrections observed in inelastic light scattering spectra by spin-density excitations 
of a slightly occupied first-excited subband. 



The existence of a spin-polarized phase of the homoge- 
neous electron gas at densities between 10 20 and 10 18 
cm and at zero magnetic field, although predicted 
years ago Q and later well established from quantum 
Monte-Carlo calculations Q, remains experimentally un- 
confirmed. It has been suggested that a partially polar- 
ized fluid state might be at the origin of the high- 
temperature weak ferromagnetism observed recently in 
La-doped calcium hexaboride U In this case, however, 
the itinerant magnetism is still a matter of debate, for 
there might be an alternative explanation in terms of 
point defects Modulation-doped semiconductor het- 
erostructures offer a unique opportunity to study the 
many-body behavior of dilute electron systems due to 
the extremely long mean-free paths that can be achieved 
with current growth techniques. Experiments carried 
out so far, using high-mobility electron gases formed 
in GaAs/AlGaAs samples, exhibit, upon depletion,the 
metal-insulator transition at about 2 x 10 10 cm -2 0,0,0 
due to the poor screening of the disorder potential of the 
ionized remote donors. In two dimensions and without 
many-body screening, this transition should take place 
for any degree of disorder |9(. Localization effects thus 
hamper further search for low-density phases of the ho- 
mogeneous two-dimensional electron gas (2DEG). Here 
we followed a different approach by raising the Fermi level 
in a single quantum well (SQW) structure in order to ob- 
tain a slightly populated excited subband. By means of 
photoluminescence (PL) spectroscopy and inelastic light 
scattering we were able to assess the properties of the 
dilute component corresponding to the second subband 
but having the 2DEG always in a metallic state because 
of the charge present in the ground-state subband. 

This Letter reports on possible ferromagnetic order in 
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a modulation-doped GaAs SQW, as predicted by self- 
consistent density-functional calculations using an exact 
exchange potential for the 2D electron system 01 an d 
as determined from PL and Raman measurements at low 
temperatures and excitation powers. Exchange terms of 
the Coulomb interaction lead to a thermodynamical in- 
stability of the 2DEG formed in the quantum well, as 
soon as the Fermi energy becomes resonant with the bot- 
tom of the first-excited subband 11]. At this first-order 
phase transition the electron gas component correspond- 
ing to the excited level breaks into spin-polarized do- 
mains with different in-plane magnetization, as expected 
for planar ferromagnets to minimise the stray field. Ev- 
idence for parallel spin alignment is obtained from PL 
measurements of the light polarization at low magnetic 
fields and from inelastic light scattering by spin-density 
excitations involving electronic transitions between the 
first and second excited subbands. 

A 2DEG with a mobility of 8 x 10 5 cm 2 /Vs and a den- 
sity of 6 x 10 11 cm -2 at 4.2 K forms in a modulation- 
doped GaAs/Alo.33Gao.67As quantum well grown by 
molecular-beam epitaxy. We present results of two sam- 
ples with 25 and 30 nm well thickness. The 2DEG is 
contacted from the surface by In alloying in order to ap- 
ply a dc bias between it and a metallic back contact. For 
the 25-nm-wide well and without bias only the lowest 
subband is occupied with Fermi energy Ep rs 25 meV. 
The energy separation to the second subband is Eq% w 28 
meV. For the wider well sample there is a small occupa- 
tion of the first-excited level at zero voltage. Photolumi- 
nescence spectra were excited using low power densities 
in the range of 0.5 to 5 W/cm 2 and recorded with opti- 
cal multichannel detection. Magneto-PL measurements 
were carried out in Faraday as well as Voigt geometry by 
placing the sample into the cold bore of a 5 T split-coil 
magnet. Inelastic light scattering spectra were recorded 
in backscattering geometry using a line focus in order to 
avoid heating of the 2DEG at the laser powers required. 
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Band-structure calculations of the dilute 2DEG phases 
were performed in the framework of density-functional 
theory (DFT) ^1 ■ Implementation of DFT requires non- 
trivial approximations in the exchange and correlation 
contributions to the total energy. In the widely used 
local-density approximation (LDA), both exchange and 
correlation potentials correspond to results obtained for 
the homogeneous 3D electron gas. A key ingredient in 
our calculations is the exact treatment of the exchange 
energy whereas the remaining (small) correlation 

contribution to the total energy is calculated within LDA. 
We also approximately accounted for low-temperature ef- 
fects in the exchange potential. Results obtained by solv- 
ing the self-consistent Kohn-Sham equations correspond- 
ing to the experimental sample geometry are summarized 
in Fig. 1. For a given electric field which mimics the ap- 
plied bias, we calculate the total density n and a set of 
energy eigenvalues Ef a and eigenfunctions £f Q , where 
i=0,l is the subband index, a=PI, PII, F the phase in- 
dex and cr =t, I the spin component. Depending on total 
density there are three possible configurations: a) The 
paramagnetic phase denoted PI, where only the ground 
subband is occupied with equal population for spin-up 
and down components. This configuration is stable until 
E^ 1 = Eqx,pi = E x .pi - E 0> pi. b) Phase PII, with both 
the first and second subbands paramagnetically popu- 
lated. This configuration is a self-consistent solution at 
high densities, where Ep 11 = E i,pn/2 + Trh 2 n/2m* , 
where m* is the electron effective mass, c) A ferromag- 
netic phase F for which the lowest subband remains para- 
magnetic but only one spin component of the second sub- 
band (say f ) is occupied and the other empty. This solu- 
tion exists for Ep > E^ F but collapses towards the PII 



configuration when E F = EL p because of the kinetic 
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energy cost to maintain a single spin population in the 
second subband. 



Figure 1 shows two remarkable effects, both beyond the 
reach of LDA: The abrupt decrease of the intersubband 
spacing Eoi at about 0.72 x 10 12 cm -2 when transforming 
from phase PI to PII and the very existence of the spin- 
polarized phase F. The sudden renormalization of Eqi, 
which is a consequence of intersubband exchange terms, 
has been discussed in Ref. For an incipient popula- 
tion of the excited subband the leading exchange contri- 
bution to the 2DEG compressibility is negative, such that 
the electron system becomes thermodynamically unsta- 
ble against the shift of a macroscopic amount of charge 
from the ground to the first-excited subband. When 
Ep = Eqi electrons start to be transferred from states 
at the Fermi surface of the lowest subband to Brillouin- 
zone center states in the excited 2D level. This leads 
to the discontinuities observed in PL associated with the 
spontaneous breaking of translational invariance. 

In contrast, spin polarization is favored by intrasub- 
band exchange terms at low densities of the excited sub- 
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FIG. 1: Dependence on total density n of the intersubband 
energy _Eoi calculated within exact-exchange theory for a 25- 
nm-wide SQW at 5 K. Different symbols correspond to the 
paramagnetic solutions with one (PI) and two (PII) occu- 
pied subbands and the ferromagnetic phase (F), as indicated. 
Lines indicate the position of the Fermi level in each phase. 
The inset shows a sketch of the conduction band profile of the 
SQW. 



band n± = nil + riif, yielding a negative non-linear con- 



tribution that goes as 
it is given by 
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Thus, the system gains energy by unbalancing the spin 
population of the second subband. The last term in Eq. 
(1) is an exchange- hole energy gain in the excited sub- 
band, which reduces the overall Coulomb repulsion in the 
2DEG due to the Pauli principle. With increasing m, 
positive contributions from the kinetic and Hartree ener- 
gies overcome the intrasubband exchange part and the F 
phase collapses towards the paramagnetic PII phase. 

Direct evidence of the existence of the F phase is ob- 
tained from PL spectra recorded at low temperatures 
using very low laser excitation. Figure 2 shows a set 
of spectra of the 25-nm-wide SQW taken at 1.7 K and 
3.5 /iW laser power with different gate voltage. As the 
Fermi level is raised towards resonance with the second 
subband, an additional peak becomes apparent at about 
1.2 meV below the energy E\ of optical recombination 
processes from the first excited electron subband. The 
F feature literally pops up in PL, increasing in inten- 
sity while merging together with the E± peak before its 
abrupt renormalization takes place. We emphasize that 
peak F is only observable if the lattice temperature is 
kept below 5 K and/or the power density does not ex- 
ceed 10 W/cm 2 . On the contrary, at higher excitation 
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FIG. 2: PL spectra of the 25-nm-wide SQW in the energy 
region of the E\ transition for different gate bias V 9 taken 
at 1.7 K with 0.5 W/cm 2 power density. As an example, the 
inset displays a full-range PL spectrum (solid curve) and the 
fitting function (dashed line) used for lineshape analysis 0. 



levels the E\ peak instead of disappearing displays just a 
slight line broadening and a strong increase in intensity. 
We also note that close to the PI-PII phase transition 
the 2DEG is very unstable, a situation which is reflected 
in the uncontrolled change of the Fermi level, and thus 
of the PL spectrum, even without varying the external 
voltage. In fact, the spectra of Fig. 2 correspond to 
snapshots of 1 s during cycles where the electron system 
goes back and forth having successively different occu- 
pations of the second subband. The first order of the 
transition manifests itself in the transformation of the 
Ei lineshape from Lorentzian to Gaussian and back to 
Lorentzian during the process of renormalization. Such 
a change from homogeneous (lifetime) line broadening to 
an inhomogeneous one (level distribution) is consistent 
with the formation of density domains in the second sub- 
band during the phase transition. The F feature is then 
associated with the recombination of electrons within the 
domains, whereas the peak E\ arises from hot lumines- 
cence coming from regions of the SQW plane where the 
second subband is still empty. 

The energy position of the E\ and F peaks measured 
at 2 K with 20 /iW is plotted in Fig. 3 as a function 
of the Fermi energy. When the Fermi energy equals E±, 
an abrupt subband renormalization by about 2 meV oc- 
curs. Afterwards, one finds the 2DEG in the paramag- 
netic phase PII and the variation of the subband energy 
with bias, i.e. Fermi level, is continuous but exhibits a 
clear hysteresis. The stability range of the F phase com- 
prises Fermi energies only within 3 meV below E\, and 
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FIG. 3: Energy position of the Ei and F peaks referred to 
Eo, the bottom of the ground subband, versus Fermi energy. 
The corresponding electron temperature is 10 K. Solid (open) 
symbols correspond to data taken with increasing (decreasing) 
bias. Solid lines are a guide to the eye. The inset shows the 
phase diagram of the 2DEG at zero magnetic field obtained 
from PL measurements. Grey data points are results of Ref. 
p"H obtained at high power densities. 



in this region the energy of the F peak is almost inde- 
pendent of the Fermi level. This is a very important 
point, for it rules out an explanation of our observations 
in terms of the Fermi-edge singularity [13j . 

The inset to Fig. 3 illustrates the phase diagram of 
the 2DEG conceived from the sets of data obtained at 
different temperatures and laser powers. Effective elec- 
tron temperatures were determined from the width of the 
Fermi edge in PL spectra. Below a sample temperature 
of 25 K, the temperature of the 2DEG is mainly given 
by the excitation-power level. The lowest value of 5 K is 
obtained with the sample immersed in superfluid helium 
and using only around 1 fxW of laser power. The region 
in the phase diagram representing the stability range of 
the F phase corresponds to the Fermi levels and temper- 
atures for which the peak F is apparent in PL. At 13 K, 
where the F phase collapses into the PII one, there is a 
triple point. We note that this temperature corresponds 
to the formation energy of the F phase, as given by the 
difference between E\ and F peak energy. As reported 
earlier [11|, there is a critical point at Tc=35 K. For com- 
pleteness we point out that the same phenomenology has 
been observed in three other SQW samples with different 
well widths. 

Having experimentally proven the existence of an ad- 
ditional phase of the 2DEG, we now present the evidence 
for ferromagnetic order stemming from magneto-PL mea- 
surements. We consider first the case of magnetic fields 
applied perpendicular to the plane of the SQW (Fara- 
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day geometry). Whereas the E\ peak splits into two 
Zeeman components with clearly different circular po- 
larization er+, cr~ with the expected gyromagnetic fac- 
tor [l4|, the F peak does not show any splitting at all 
nor any degree of circular polarization. On the contrary, 
its intensity decreases linearly with field, disappearing at 
about 1.5 T. Such behavior is compatible with the as- 
sumption that the magnetization within the density do- 
mains is randomly, but in-plane, oriented. A magnetic 
field in the perpendicular direction destroys the domains 
since Bj^ significantly reduces the wave-function overlap 
in the second subband and thus the intrasubband ex- 
change interaction. 

For magnetic fields applied along the 2DEG plane the 
intensity of the F peak remains essentially constant. In 
contrast, its degree of linear polarization defined as 7 = 
(Iv — Ih)/{Iv + Ih) increases linearly with B|| by 30% 
in the range from zero to 0.6 T. Here ly and Ih are the 
intensities of the F peak for vertical and horizontal linear 
polarization, respectively, normalized to the intensity of 
the PL spectrum in the energy region just underneath 
the Fermi edge, where the emission is expected to be 
independent of its polarization. The field is aligned along 
the horizontal direction, but rotated 90° with respect to 
the direction of light collection (Voigt configuration) . In 
this geometry, light emitted with circular polarization 
along the field is detected as being vertically polarized. 
Hence, the observed linear increase in 7 can be ascribed 
to the successive orientation of the domain magnetization 
parallel to the external field. Thus, 7 is given by the 
Brillouin function Bj which in the low-field limit yields 
0] 7 = (>/+ l)ffA t s-B||/(3fc_BT e //). Here is J the average 
total spin of the domains in units of h/2, g =-0.44 is 
the gyromagnetic factor of GaAs 0], \ib is the Bohr 
magneton, ks is the Boltzmann constant, and T e // is 
the electron temperature. From the measured slope for 
7 we infer an average number of 1 /2-spins, i.e. electrons, 
within a domain of about 50 (25 nm SQW). For a typical 
electron density of 5 x 10 10 cm -2 corresponding to the 
jump at the PI-PII phase transition the mean radius of 
a domain should be around 200 nm. 

Another striking result which speaks for a spin po- 
larization of the dilute electron system in subband '1' 
is obtained from inelastic light scattering. Collective 
spin-density (SD) excitations of the 2DEG are known to 
be shifted down in energy from the intersubband single- 
particle (SP) gap by many-body effects, the so-called ex- 
citonic shift. It represents the vertex correction due to 
exchange-correlation terms of the Coulomb interaction 
The determination of its magnitude is straightfor- 
ward from the energies of the excitations Tiujsp and husD 

measured in depolarized light scattering spectra, accord- 
■ — 1 2 _ 2 

ing to Hi : -2-y S Dn 2 D = " 9 " S p SF ■ Here -y S D represents 
the exchange-correlation energy and n<iD is the electron 
density of the occupied subband. Figure0|shows the val- 
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FIG. 4: Exchange- vertex correction to the spin-density excita- 
tion involving electronic transitions between subbands '1' and 
'2' as a function of the intersubband spacing between the two 
lowest subbands of a 30-nm-wide quantum well. Open(closed) 
symbols correspond to measurements using low(high) power 
densities, i.e. the effective electron temperature being 5 and 
15 K, respectively. Curves are a guide to the eye. 



ues of the exchange vertex correction to the spin-density 
excitation associated with electronic transitions between 
the first and second excited subbands of a 30-nm-wide 
SQW. In the graph, the first-excited subband is being 
filled for decreasing subband spacing E i- Open symbols 
correspond to data obtained at a very low power density 
such that the F peak is apparent in PL. For comparison 
the full circles represent the high-power case, where the 
effective electron temperature is too high for the spin- 
polarized phase to exist (peak F is no longer observ- 
able). The exchange vertex correction exhibits a strong 
enhancement at low electron temperatures and very low 
filling of subband T',i.e. in the stability range of phase 
F, whereas at high power levels it extrapolates to almost 
zero. The results are independent of laser power, when 
the renormalization of E\ sets in. As demonstrated ear- 
lier for fractional quantum Hall states with odd filling 
factor [l?! and in GaAs quantum wires in magnetic fields 
(lif , this is clear indication of spin polarization. The cost 
in exchange energy to flip a spin in a ferromagnetically 
ordered state exceeds many times that of the nonmag- 
netic case. We emphasize that such an enhancement of 
vertex corrections is clearly observed for the 1-2 elec- 
tronic transitions rather than the '0, 1' subbands, for in 
the latter case most of the electrons involved in the spin 
excitation are in a paramagnetic state. 
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